Simulation of many-body interactions by conditional geometric phases 
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It is shown how to exactly simulate many-body interactions and multi-qubit gates by coupling finite 
dimensional systems, e.g., qubits with a continuous variable. Cyclic evolution in the phase space of 
such a variable gives rise to a geometric phase, depending on a product of commuting operators. The 
latter allows to simulate many-body Hamiltonians and nonlinear Hamiltonians, and to implement a 
big variety of multi-qubit quantum gates on both qubits and encoded qubits. An application to the 
quantum amplitude amplification algorithm will be discussed. 
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I. INTRODUCTION 

One of the most challenging and interesting prob- 
lems in theoretical physics concerns the understanding of 
quantum many-body systems. In view of the exponen- 
tially many degrees of freedom involved, it is generally 
agreed that simulation of such systems is a task that can- 
not be efficiently tackled by any classical computers. On 
the other hand, as suggested by Feynman Q the growth 
in computational resources is only linear on a quantum 
computer which is itself a quantum many-body sys- 
tem. Therefore even a modest size quantum computer, 
e.g., containing only a few tens of qubits, could outper- 
form a classical computer. We should like to stress that 
the feasibility of such a "devoted" quantum computer 
problem is expected to be definitely greater than the one 
of a general purpose quantum computer [QJ|]. 

Simulation of Hamiltonians has very practical applica- 
tions in quantum control [|J and quantum information 
theory ||. A typical goal is to achieve, i.e., simulate, as 
many as possible desired entangling Hamiltonians start- 
ing from a given one and the ability to perform local 
operations []t3| — |l Of] or more generally to enact transfor- 
mations drawn from a given "natural" set of available 
interactions. 

If we have two Hamiltonians Hi and H2 we can sim- 
ulate the Hamiltonian X\Hi + X2H2 (Ai and A2 are real 
numbers) and the Hamiltonian i[H\,H2\ | pT| . Thus any 
Hamiltonians in the algebra generated by Hi and H2 can 
be simulated. This simulation method does not require 
that the two Hamitonians Hi and H2 commute but re- 
quire that one can alternately turn on them for arbitrar- 
ily short time over infinite steps. Physically switching 
Hamiltonians over arbitrarily short time is not a easy 
task and unrealistic in some sense. We shall propose a 
method to do exact Hamiltonian simulations over finite 
steps. 

In this paper we consider a "hybrid" physical sys- 
tem consisting of N qubits coupled with a continu- 
ous variable system. The state space is then given by 
U = (<C 2 )® N <g> /ioo, where h m := span{|ra)}£L , is the 
standard Fock space for a single mode described by the 
bosonic annihilation and creation operators a and , re- 



spectively. We assume the following Hamiltonians Jl2 13 
to be realizable 



Hi = —iX(aa f — a* a) A 
H2 = uja) a A, 



(1) 
(2) 



The operator A is a pure qubit operator that in the sequel 
will be mostly either a Pauli operator cr ja (a = x, y, z) 
for ion i or a collective angular momentum operator 
J a = \ J2iLi a ia- The case of A = J a corresponds to N 
ions with each driven by identical Raman lasers. From 
the Hamiltonians above we can have the conditional dis- 
placement operator and the conditional rotation operator 



Ui=e 
U 2 



— iH\t — XtA(aa^ — a* a) 



-iH 2 t 



(3) 
(4) 



respectively. These two operators will play a crucial role 
in this paper in which we will try to address the follow- 
ing question: given these two kinds of physical operators, 
what is the maximal set of Hamiltonians that we can sim- 
ulate? In Ref. jljj it has been shown that a conditional 
displacement of the vibrational mode of trapped ions can 
be used to simulate nonlinear collective and interacting 
spin systems including nonlinear tops and a universal 
two-qubit gate, independent of the vibrational state of 
the ion. The scheme in JTJ] has been further extended in 
order to realize the nonlinear Hamiltonian J 2 [l5| ] and 
multi-qubit quantum gates fl6|| . 

A general framework for quantum information process- 
ing (QIP) with hybrid quantum systems has been pro- 
posed in Ref. [jl7|. There it has been proved that an 
universal set of Hamiltonians for hybrid quantum compu- 
tation is provided by the following: {±a x x, ±a z x, ±a z p}. 
The position and momentum operators x and p satisfy 
the canonical commutation relation [x, p] = i. The 
ability to turn on and off the Hamiltonians from this 
set allows one to produce conditional displacements that 
in turns allow to enact Hamiltonians that are arbitrary 
polynomials of the a Xl a y ,a z ,x, and p. We can use an al- 
ternative set of Hamiltonians {a x a'a, a z a^a, x}. This set 
can be easily reduced to the above one as follows: 
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± a x x = e^ a *^ a p e** (5) 
±a z x = e^°^ a p e ±f ^ ata (6) 
±a z p = e ±] i' y ^ a x e^ a ^ a (7) 

The operator p can be obtained from (^) when the spin 
is in the state |0) or |1). These equations show that not 
only conditional displacements but also conditional ro- 
tations are useful in quantum computers. We will use 
both the conditional displacement operators and condi- 
tional rotation operators to simulate many-body interac- 
tion Hamiltonians and implement various quantum gates 
@|| and algorithms 

Our scheme is interestingly related to one kind of ge- 
ometric phase, i.e, the geometric phase in phase space 
plLEgl . However now our geometric phase is depen- 
dent on some operators, i.e., it is a conditional geomet- 
ric phase. Recently both the adiabatic |^3| and non- 
adiabatic geometric phases have been suggested as 
potential candidates for realizing quantum computers 
displaying some built-in fault tolerant features pB| p9[ . 
We shall see that even geometric phases in phase space 
could be useful for QIP. 

The paper is organized as follows. In Sec. II we first 
review the concept of geometric phase in phase space. 
In Sec. Ill we show how to generalize this phase to 
'conditional' geometric phase and simulate many-body 
interactions. Once these Hamiltonians are available, the 
implementation of some quantum gates become somehow 
straightforward. We discuss the implementation of quan- 
tum gates on both qubits and encoded qubits in Sec. IV. 
In Sec. V we first show how to realize general projec- 
tors and then how to implement the quantum amplitude 
amplification algorithm using the projectors. The con- 
clusions are given in Sec. VI. 

II. GEOMETRIC PHASES IN PHASE SPACE 



D(0) = 1, and V(-a) = P(a) t . From Eq (g) it follows 
also the useful identity 

eJ im03 Q *) = T>{-p)V{-a)V{p)V(a). (10) 

We consider now a closed loop 7 which is an N— 
sided polygon with sides ay (j — 1,2,.. .,7V) such that 
^2^=1 a j = 0- The total transformation associated with 
7 is given by Q 

D 7 = V(a N ) ■ • • D(ai). (11) 

An arbitrary closed loop 7 can be approached in the 
limit of N — > 00. In this case the total transformation is 
just a phase factor 2? 7 = exp(iO), where 

= i <j){xdp~ pdx). (12) 

The phase O neither depends on the form of the loop nor 
on the speed of the transformation but just on the area 
of the loop. For this reason it deserves the name of a 
geometric phase. 

III. SIMULATION OF MANY-BODY 
INTERACTIONS 

Let us recall the Hamiltonian simulation ||. There is 
a set of Hamiltonians {Hi} and a class of allowed opera- 
tions like local unitaries or local operations and classical 
communications. The aim is to produce a desired evolu- 
tion e~ lHt from the set of the Hamiltonians and the class 
of allowed operations, where H is the simulated Hamil- 
tonian the t is the simulated time. 

Sometimes the class of local operations is not cheap 
resource. In this paper we will not use them and just 
begin from a set of Hamiltonians to do the simulation. 



We begin by a brief review of the geometric phase 
within the formalism recently given by Lius ]22] | . Phase- 
space translations are represented by the displacement 
operator 

V (a) = D (V2aj = e*^ - ** , (8) 

where the complex quantity a — x + ip parametrizes the 
displacement and D (a) = exp(aa^ — a*a) is the usual 
displacement operator in quantum optics. These opera- 
tors satisfy the relation 

V(P)V(a) = ei lia{fla * ] V{a + /?), (9) 

If a = xx + ipi and f3 — x-i + ip2 one has Im(/3a*) = 
xip 2 - x 2 pi- 

Eq.(||) tells us that the displacement operators real- 
ize an unitary (projective) representation of the addi- 
tive group of complex numbers. In particular one has 



A. Nonlinear Hamiltonians and two— body 
interactions 

A conditional displacement operator reads T>(Aa), 
where the Hermitian operator A usually represents a dis- 
crete variable observable. By using the spectrum decom- 
position of the operator A, A — Y2k=i ^k\k)(k\ the dis- 
placement operator can be written as 

N 

V(Aa) = ^2\k)(k\V(X k a). (13) 

k=l 

From the above equation one can see clearly that the 
amount of displacement is dependent on the eigenvalues 
of A. By replacing the displacement operators in Eq. (]To|) 
with conditional displacement operators, we immediately 
obtain 
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J8AB 



= V(-Bf3)V(-Aa)V{Bf3)V(Aa) , 



(14) 



where 8 =Im(/3a*) and operator B commutes with A. 
Now the geometric phase depends on the operator prod- 
uct AB. Therefore we can simulate the Hamiltonian AB 
if operators A and B commute and belong to one system. 
For example if A = B = J z , we get the nonlinear Hamil- 
tonian J 2 , which is useful for creating GHZ multipar- 
tite entangled states |3C]]3l[] and generate spin-squeezed 
states |32j . If A and B belong to different systems we may 
simulate the two-body interaction Hamiltonian A ® B. 



B. Three— body interactions 

In the context of NMR quantum computation, Tseng 
et al. 1 33 1 proposed a method to simulate three-body 
interaction of the type <j z ®o z ®<j z . Due to the self-inverse 
property of a z simulation of the three-body Hamiltonian 
becomes relatively easy by two-body interactions. The 
simulation can be obtained in the following way 



gi^CTla-lgffSs gif<Tl !C ®a'2a 



(15) 



Now we show how to create 'nonphysical' three-body 
interaction Hamiltonians H — A ® B ® C for general 
operators A, B, and C which act on three different sub- 
systems, respectively. To this aim let us introduce the 
operator 71(8) — exp(i6*a^a), it entails a rotation in the 
phase space 



TZ(-9) 



P 



11(0) = 



COS ( 

sin t 



— sm( 
cos 8 



Directly from the above equation we obtain 
V(ae w ) =K{0)V(a)n{-d). 



(16) 



(17) 



The conditional rotation operator is given by 1Z(8C), 
which makes a rotation in phase space conditioned on 
the operator C. As in Eq.(|17|) we get 



(18) 



V{ae i9C ) = TZ(8C)V(a)TZ(-8C). 
Then we can replace a with ae t9C in Eq.(14) and obtain 

e -iTA®B®sin(6C+rf>) 

= V(-B[3)V(-Aae i96 )V(Bf3)V(Aae wd ), (19) 

where r = \a/3\ and <f> = ar g(«) — ar g(/3). Eq.([l9|) tells 
us that we can simulate the three-body Hamiltonian 



H = \A®B® sm(8C + <j>) 

by the following sequence 

V(-BP)K(6C)V(-Aa)1l[-6C) 
xV{Bf3)TZ(8C)V(Aa)TZ(-8C). 



(20) 



(21) 



Therefore the 'nonphysical' three-body interaction is ex- 
actly achieved by eight two-body physical interactions. 
Specifically we obtain the Hamiltonian 



H = XA <g> B <g> sin(0C). 



(22) 



If we make a small rotation, i.e., the angle 8 is small 
enough, we can approximately realize the three-body in- 
teraction Hamiltonian 



H = X'A ® B <g> C. 



(23) 



among three operators A, B and C, where A' = A#. Re- 
markably for some operator C we can even exactly realize 
the three-body interaction. 

Let us consider an operator C which satisfies either 
C 2 = 1 (self-inverse operators) or C 2 = C (idempo- 
tent operators). For such an operator we have sin(#C) = 
sin 8C. Therefore Eq.(|2|) shows that we can exactly real- 
ize the Hamiltonian A®B(g)C. The self-inverse operators 
we often encounter are Pauli operators and parity oper- 
ators et al. Most quantum gates in quantum computer 
are also self-inverse, such as Hadmard gate, controlled- 
NOT gate (C NOT ) Q, SWAP gate and Toffoli gate 
Examples of idempotent operators are of course given by 
projection operators. 

Summarizing the three-body interaction Hamiltonian 
A eg) B ® C is approximately simulated for any opera- 
tors A, B, and C, and is exactly simulated if one of the 
operators is self-inverse or idempotent. For the case of 
A = B = C = J Zl we approximately have the Hamilto- 
nian J z (g) J z (g) J z . By choosing the self-inverse C = a z , 
we exactly have the Hamiltonians J z <B> J z <E> c z , and 
a z ® o z <g> a z . 

Both Eq.([l4]) and Eq.(^l|) correspond to a close loop 
(a parallelogram) in phase space. Now the resulting geo- 
metric phases are dependent on a product of commuting 
operators, which are just the desired simulated Hamilto- 
nians. 



C. Many— body interactions of qubits 

In this subsection we address the simulation problem 
of some many-body interactions of qubits which are re- 
lated to quantum spin models in condensed matter the- 
ory. From Eq.(20) we have the Hamiltonian 



Let 8 



H c = Acos(#J Q ), 
H s = Xsm(8J a ). 

= 7r we obtain 

H c = Xcos(nN/2)a a ® a a < 
H s = Asin(7riV/2)(T Q ® a a ( 



) (T Q , 



(24) 
(25) 



(26) 
(27) 



From Eq.(||) with even N and Eq.@ with odd N, we 
can simulate, for arbitrary N, the multi-qubit Hamilto- 
nians 
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H 



«± 



±A(T n 



(28) 



The above many-body Hamiltonians are very useful. For 
instance, let the initial state the ./V-qubit system be 
lO) 8 ^ and evolve according to the Hamiltonian H x+ , 
then the state vector at time t is easily obtained as 



9,N) = cos(0)|O) 



isin(0)|l 



(29) 



where 9 = Xt. The state 1 7r/4, N) is a multipartite entan- 
gled state, which is known to play an important role in 
quantum information theory. 

From Eq.(p8|), we have the Hamiltonian <J\ a oia(pi = 
x,y, z) for N = 2, which are just the Ising interactions. 
As (J\ x 02x, viyV2y, and o\ z 02z commutes with each other, 
we may simulate a general two-qubit Heisenberg model 

M 



Even for noncommuting set { XiHi}, we still have a chance 
that a decomposition similar to Eq.(|36|) exists (37|. For 
instance, the operators in Eq. (|33|) , a a — a® N ', act as the 
encoded Pauli matrices. Then we have the identity 

exp(— i(f>[cos Ba z + sin9a y )/2] 
= exp(i6a x /2) exp(—i(j)a z /2) exp(—i6a x /2), (37) 



which implies that we can simulate the Hamiltonian 
H = X z a z + X y a v = X z af N + X y af N . 



(38) 



for odd N = Am + 1 using a x = af N and b z = af N . 



IV. SIMULATION OF QUANTUM GATES 



H — X x U\ x Ui x + Xy<Ji y G2y + X z a\ z (J2z 



(30) 



For N > 2, i.e., the many-body case, we already have 
the Hamiltonians erf N , af N , and of N , which satisfy the 
commutation relation 



JgliV 



y J 



•JVf 



1±(-1)> 



(31) 



Here the subscripts + and — indicates the anticommu- 
tation and commutation, respectively. From the above 
equation we see that the three operators either commute 
or anticommute. For even TV, any two of the three oper- 
ators commute with each other. Then we can simulate 
the following Hamiltonian |3l| 



H = X x a®W + X y af 2k ^ + X z af^ 
(k = 1,2,3...). 



(32) 



However for odd N, any two of the three operators 
anticommute but do not commute with each other. For 
N = 4m + 1 (m = 0, 1, 2, ...), from Eq.@, we obtain 



\af N ,ai 



= 2ia? N . 



(33) 



Thus the operators af N /2, af N /2 and af N /2 realize the 
su(2) Lie algebra. For odd N = Am + 3 (m = 0, 1, 2, ...), 
from Eq. (|3l|) , we obtain 



y 



(34) 



In this case of N /2, af N /2 and -af N /2 realize the su(2) 
Lie algebra. For even iV we can find three operators 



r ®(2fc) 



5(2fc— 1) 



/2, a z ® I 



<8(2fc-l) 



12 (35) 



satisfy the su(2) commutation relations. These realiza- 
tions of su(2) Lie algebras will be used in the later dis- 
cussions of quantum gates on encoded qubits. 

From any set of commuting Hamiltonians {XiHi} we 
can simulate the Hamiltonian XiHi since 



cxp(^ XiHi) = ]^[exp(A 4 ffi 



(36) 



So far we have showed how to obtain, by cyclic condi- 
tional evolutions in phase space, an operator-dependent 
geometric phase and how to use these operators in order 
to simulate two-body and many-body interaction Hamil- 
tonians. Now we make use of the above-developed for- 
malism to explicitly construct some important quantum 
logic gates. We emphasize that the continuous e.g., vi- 
brational degree of freedom is only required during gating 
and acts like a databus. 



A. Two-qubit gates 

Controlled-NOT gate |Q and controlled phase gate 
(C P ) @: Let A = (1- a u )/2, B = (l-a lx )/2, and 6 = 
tt in Eq.(pl|), the controlled-NOT gate is immediately 
obtained as 



Cnot = cxp 



-i-(l - (T lz ){\ - a 2x ) 



(39) 



The first bit is the control bit and the second is the target 
bit. Similarly the controlled-phase gate is obtained as 



Cp = exp 



- a lz )(l - a 2z ) 



(40) 



A simple relation exists between the controlled-NOT 
gate and the controlled-phase gate, 



Cnot = exp 



.7T 
>—( 

4 



Cp exp 



7T 



(41) 



i.e., they differ only by local operations. 

SWAP gate : As we can simulate the Heisenberg 
Hamiltonian (|30|), the SWAP gate is easily constructed 
as 



Gswap = e ~ %1 i { ~ c ' 1:cC ' 2:c+ ' yiy ' y2y+ ' yiz ' y2z ~ 1S) 
1 

= O (1 +&lxV2x + V\y&2y + 0\ z O\ z 



(42) 
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B. Three-qubit gates 

Toffoli gate (T) @: Let A = (1 - a lz )/2,B = (1 - 
cr 2z ) /2, and C = (1 - cr 3a; )/2, in Eq.(§f), the three-bit 
Toffoli gate is obtained as 



exp 



- <T!,)(1 - cr 2z )(l - cr 3 iO 



(43) 



which is also called (controlled) 2 -NOT (C^ t) S a ^ e an( i 
important for universal computation. Note the Toffoli 
gate is exactly realized since the operator C = (1 — a^ x ) / 2 
is a projector. Such a gate can of course be also im- 
plemented by constructing appropriate networks of one 
qubit and two qubit gates Jll] . The construction showed 
here is somehow more direct. Notice also that an alterna- 
tive way to obtain the Toffoli gate is recently discussed 
in Ref. |l6|]. The implementation of the more general 
(controlled) JV -NOT will be discussed in Section V. 

Fredkin gate (F): Another example of a three qubit 
gate relevant for QIP is provided by the controlled- 
SWAP Fredkin gate 139]. The following commuting 
Hamiltonians (1 — <Jiz) o~2a o~Za{<^ = x,y,z) can 
be realized exactly (|2j). Therefore we can simulate the 
unitary operator 

which is just the Fredkin gate. 



C. Quantum gates on encoded qubits 

In order to perform universal quantum computations it 
is sufficient to be able to make arbitrary single qubit ro- 
tations together with controlled-phase gate. For encoded 
qubits one problem is how to make logical operations pCj ] 
of them and the above two logical operations on encoded 
qubits are needed. We will discuss two typical codes: the 
active error correction codes Ell and the passive codes 
on decoherence-free subspaces [|42|. 

First we consider the error correction codes with odd 
number, the linear codes proposed by Steane ifjlf . He 
have devised two encoding, the first of which protects 
only against decoherence 

|o c > = i(|iii> + |ioo) + |010) + |001)), 



Uc) = Ox 



o- x \0c 



(45) 



The encoded qubits span a representation space of 
su(2) Lie algebra generated by the operators af N ,a® N 
and -of N (§§). Here N is either 3 or 7. That is to say, 
these three operators acts on encoded qubits as encoded 
o x , <jy and a z . Then we consider the code mapping 1 
qubit into 5 qubits presented in Ref. B3[ 



|0 C ) = Q|00000), |l c ) = Qllllll), 



-(<r a 

-{<7z 



i 

Or 



<j x 1 1 1 
01® t) cyc 



1) 



eye 



cycl , 



(47) 



where the subscript "eye" indicates that addition of all 
five cyclic shifts. Obviously the commutators Jct® 5 , Q] — 
(a = x,y,z) and then erf 5 a® 5 , and crf 5 (|33|) act as 
encoded o~ Xl a y and er 2 , respectively . 

Now we consider even-number codes in decoherence- 
free subspace p2j. We can have a code of two qubits 



|0 o > = |01>, |1 C ) 
and a code of four qubits 



|10). 



|0 C ) = 10011) 
|lc) = |1100) 



|0110) + |0101), 
|1001)+ |1010), 



(48) 



(49) 



For these two cases, it is easy to see that the two- 
dimensional subspace span{|0c), |lc)} is decoherence- 
free subspace. The corresponding encoded Pauli matrices 
are given by Eq.p5|). 

The controlled-phase gate for encoded qubits is easily 
constructed as 



Cp = exp ' 

Explicitly for three-bit code we give the gate as 

- , IT 



(50) 



exp[ 



+ 1 Q 



S03 



l^+crf 6 ) 



(51) 



which can be realized by our simulation method. Then 
we have both the rotations of single encoded qubit and 
gate Cp for two encoded qubits, which are enough for 
the quantum computation on the encoded qubits. 



The second is capable of decoding with general 1-bit er- 



V. APPLICATIONS TO QUANTUM 
ALGORITHMS AND DISCUSSIONS 



71' 



11111111) + |0101010) + 11001100) 



+ 10011001) + 11110000) + |0100101) 
+ |1000011) + 10010110)), 

cr x cr x (g) a x ® a x (g> er x <j x ® o x |0 C ) 



(46) 



In this section we show some applications of the simu- 
lation schemes discussed so far to the implementation of 
quantum algorithms. 
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A. Realization of projectors 

Let us consider the Hamiltonian 

H k = Xcos[e k (J z + (j) k )}. (52) 

Since the Hamiltonians H k commute with each other for 
two different k, we have the sum of them 



(53) 



k=0 



This Hamiltonian can be obtained by N + 1 cyclic evo- 
lutions in phase space. If we choose 6 k = jfPi , A = 
-^qrj-and <pk = N/2 + n (k — 0...N), the Hamiltonian 
becomes 



H = 



N 



1 N 

+ 1 ^ 



2-Kk 
N+ 1 



(AT - n) 



(54) 



where Af = J z + N/2 whose eigenvalues range from to 
N. From the above equation and using the identity 



N 



1 N 

+ 1 ^ 

k=0 



2-Kk 
N + 1 



(n — n') 



we obtain the Hamiltonian 

H = 5j\f n 



(55) 



(56) 



which implies that we have realized the projector P n 
which project the state to the symmetric subspace with 
excitation n. For instance, for n — or N, the projectors 
are just 

P = |00...0)(00...0|, (57) 
P N = (58) 

A general Hamiltonian F(J Z ) can be written as 

N 

F{J z ) = Y J F(n-N/2)P n . (59) 

Since all P n commute we can simulate the general Hamil- 
tonian F(J Z ). 

Once we have realized the projector Pq we can realize 
any one of 2 N projectors 



\ai,Q2, ■■■,a n )(ai,a 2 , —, a n \ 



(60) 



in the iV-qubit space, where a,i £ {0, 1...N} and the state 
|eti, 02, a n ) represents that the only in the positions of 
oi, &2, and a n the qubit is in the state |1) and in the state 
|0) in other positions. The projector can be implemented 
as follows 

'Pa\a'2...a n — ® a\x@ a^x 3 ■ ■ ■ 3 ^ajyxF(j(7 aiX (T a2 x j ■ ■ ■ 1 ^a^x • (61) 



Then we realize all the projectors in the iV-qubit space. 
This can be understood from a more general point of 
view Rj. Suppose we have the Hamiltonian H and 
perform unitary operations U and XT'. Then it follows 
from the identity e - ltUHlji = Ue- ttH W that we can 
exactly simulate evolution according to the Hamilto- 
nian UHU^. Now our H and U are: H = Pq and 

U U ® a\ X&d2 X 3 '"3 &ajyx- 

The projector V ai a 2 ...a N is very useful. For instance, 
from P N := [±)(M® N , the (N + l)~bit C^ OT gate, is 
expressed as [n8| 



NOT — ex P 

= 1 — Pn + Pn cr n +ix, 



(62) 



which is a natural generalization of the controlled-NOT 
gate and Toffoli gate to many qubits. Alternatively 
one can construct the multi-qubit generalization of the 
controlled-phase gate 



U Pn = e -^ P «/ 2 



Then the C^ Q ^ gate is easily obtained as 



piJV-l 
^NOT 



Up N e l 



(63) 



(64) 



So we can straightforwardly implement C^ Q ^ gate once 
we have had the projector. 



B. Implementation of the quantum amplitude 
amplification algorithm 

In 1997, Grover presented a search algorithm J2(^] that 
identifies the single value xq that fulfills f(xo) = 1 for 
a function f(x) provided, e.g., by an oracle (all other 
arguments lead to vanishing values of the function). If x 
is an integer on the range between and N — 1 = 2™ — 1 , 
the search algorithm is able to find xq after on the order of 
y/N evaluations of the function. Grover's algorithm has 
been demonstrated on NMR few qubit systems p5|. In 
our previous paper |l6|| we have shown how to implement 
the Grover's algorithm. In the following we will use our 
method to implement a general quantum search, which is 
called quantum amplitude amplification algorithm pq ]. 

We can write a general quantum search operator 



Q = Q(A,x,<P,#) 



-AS%A 



-l c i9 



(65) 



which is at the heart of the quantum algorithm. Here 
A is any quantum algorithm that acts on the A-qubit 
system and 



S£ = l + (e^-l)|00...0)(00...0|, 
S* = l + (e«-l)5> fe )<T fc |. 



(66) 
(67) 
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The state \rk) is a marked state and the summation runs 
over all the marked states. Thus this quantum search is 
a multi-object search. The quantum algorithm contains 
a unitary transformation, two phase rotations and the 
marked states. When A is the Walsh-Hadmard transfor- 
mation, there is one marked state, and $ = (p = it, the 
quantum algorithm reduces to the usual Grover's search 
algorithm. 

Now we write the two phase rotation operators Sq and 
as exponential form as 



_ „Moo...o)(oo...o| 

On — C , 



(68) 
(69) 



As we have all the projectors (|6C|) and they commute 
with each other, the two rotation operators are then re- 
alized straightforwardly. Therefore the general quantum 
search algorithm can be implemented. In a recent paper 
p8[ the phase rotation operator = 1 — 2cos#e" ? | / r)(T| 
is introduced. This operator can be written as the expo- 
nential form as = e l< - 7r ~ 2,5 )' r ^ r ' and therefore we can 
realize it similarly as Sg and S*. 

A few comments on the above geometric scheme in this 
paper are now in order. From Eq.(|l9|) it can be seen that 
we still have freedom to make a generalization as 

e -iT :jk A j A h sixi{e j d j -ekGk+<t>jk) _ x>(Akake Wk ^ k ) 
xV{-A j a j e i6 i 6 i)V{A k a k e i0hd '')V{A j aje iS i 6 i). (70) 

Here Tjk = IctjCtjtl and <pjk = arg(ay)— axg(afe). So we can 
simulate the Hamiltonian like AjAkSm(9jCj — 9kCk + 
4>jk) which includes for arbitrary commuting operators 
Aj , Ak , Cj , and Ck ■ Then we further ask if we can ex- 
actly simulate the product of three arbitrary operators 
AjAkCi or more general one, the product of four oper- 
ators AjAkCiC m . The answer seems negative with our 
scheme. The reason is as follows. The general transfor- 



mation is given by Eq. (|1 1|) with a 



/A, 



Then 



if the transformation 7 is cyclic, we obtain a geometric 
phase factor given by 



'S 3 >fc r i k AjA k sin(9 j Cj-8 k Ck + tt'jk) 



(71) 



Hence we can not exactly achieve the products of three 
or more commuting operators. 



VI. CONCLUSIONS 

In this paper we have discussed a geometric scheme 
to simulate many-body interactions and to implement 
multi-qubit gates. Our strategy is based on conditional 
cyclic evolutions in the phase space of a continuous quan- 
tum variable coupled to discrete systems. The cyclic evo- 
lution leads to a conditional geometric phase factor con- 
taining the desired operators acting non trivially on the 



discrete factor. To the use of the reader let us list the 
main results of this work: 

1. We can exactly simulate two-body Hamiltonians 
like A<S>B, the three-body interaction Hamiltonians like 
A® B <g) sm(6C + 4>), A® B®C with one of them is self- 
inverse or idempotent, and many-body Hamiltonians like 
<j a ®a a ®...®<J a . Some quantum spin models in condensed 
matter theory are exactly simulated. We approximately 
simulated the Hamiltonian A ® B ® C for three operators 
A, B, and C. 2. We can simulate the nonlinear Hamil- 
tonian Jl and the more general nonlinear Hamiltonian 
F(J Z ). 3. Nearly all the quantum gates proposed until 
now, especially multi-qubit gates and quantum gates on 
encoded qubits, can be implemented with our scheme. 
4. We simulate the projectors as quantum Hamiltonians 
and implement the quantum amplitude amplification al- 
gorithm. 

In conclusion we would like to stress that, in the simu- 
lations discussed in this paper, the Hamiltonians involved 
are given only by conditional displacement operators and 
conditional rotation operators. These operators can be 
realized experimentally, e.g., in ion-traps. Therefore we 
believe that the simulation strategies discussed in this 
paper have direct practical relevance. 
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